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Abstract 



> 

' We study the non-perturbative renormalization group flow of higher-derivative gravity 

£SJ ■ employing functional renormalization group techniques. The non-perturbative contributions 

to the /3-functions shift the known perturbative ultraviolet fixed point into a non-trivial fixed 
point with three UV-attractive and one UV-repulsive eigendirections, consistent with the 
| asymptotic safety conjecture of gravity. The implication of this transition on the unitarity 

problem, typically haunting higher-derivative gravity theories, is discussed. 

X 

.«?. 

Among the many approaches to quantum gravity, a special place is occupied by higher- 
derivative gravity, which, besides the Einstein-Hilbert term, also includes fourth-order operators 
in the action. Indeed, the higher-derivative propagators soften the divergences encountered in 
the perturbative quantization, rendering the theory perturbatively renormalizable [1] and asymp- 
totically free at the one-loop level [21 El IH El E] • Unfortunately, the extra terms responsible for 
the improved UV behavior also induce massive negative norm states [7], so-called "poltergeists", 
which led to the belief that the theory is not unitary. Several arguments suggest that this short- 
coming can be cured by quantum effects [21 E], but the lack of non-perturbative methods has 
made it hard to substantiate such claims. 

Recently, the question of renormalizability has received renewed attention due to mounting 
evidence in favor of the non-perturbative renormalizability, or asymptotic safety (AS), of gravity 
[Si [TUl [TT1 [T2] . In this scenario, the ultraviolet (UV) behavior of the theory is controlled by a 



non-Gaussian fixed point (NGFP) of the renormalization group flow, with a finite number of 
UV-attractive (relevant) directions. This fixed point is supposed to provide a continuum limit 
which is predictive and safe from divergences. While initially found only for gravity in 2 + e 
space-time dimensions [13|. 114] or in the 1/N expansion |15[ I16], a candidate for such a NGFP has 
been recently observed in four dimensions thanks to the application of functional renormalization 
group equations (FRGE) to gravitational actions of the /(fl)-type [I71IIgllI51l20|l2Ill52|l551l231. 
Further support for AS comes from lattice simulations \25\ [26] . 

Since higher-derivative terms play a central role in the renormalizability and unitarity of 
gravity theories in the perturbative setting, it is of great interest to study their impact on the 
AS scenario. To this effect, in this letter we investigate the RG flow of gravity including all four- 
derivative interactions and present, for the first time, non-perturbative results in this setting. In 
particular, our work augments the previous f(R) findings by including the second power-counting 
marginal coupling related to the Weyl-squared interaction. This constitutes an important step 
in the program of studying the non-perturbative RG flow of gravity in a derivative expansion, 
since, besides higher-derivative vertices studied previously, the approximation also includes a 
four-derivative kinetic term for the graviton propagator. 

Our main result is the existence of a NGFP with three UV-attractive and one UV-repulsive 
eigendirections. Together with the condition that this fixed point provides the continuum limit for 
gravity, the latter implies a relation between the coupling constants of the theory at high energy. 
Tracing the RG flow to the IR, this translates into a relation between the four-derivative couplings 
in the effective field theory which could, in principle, be tested experimentally. Furthermore, the 
transition from asymptotic freedom to asymptotic safety may realize an old proposal [21 [8] for the 
removal of the poltergeists from the spectrum of the theory, providing the exciting perspective of 
asymptotically safe quantum gravity being unitary. 

The key tool in our investigation is the effective average action and its FRGE [30] (see, 
e.g., [281 \10\ [TT] for reviews). This setup provides a continuum analogue of Wilsons lattice 
renormalization group, realizing the idea of integrating out all fluctuation modes with momentum 
larger than a certain cutoff k (high momentum modes), and taking them into account through a 
modified dynamics for the remaining fluctuations with momentum smaller than k (low momentum 
modes). This is implemented by means of a cutoff TZ^ which, for a given k, suppresses the 
contributions of the low momentum modes, so that only the high momentum modes are integrated 
out in the path integral. The functional then defines an effective field theory valid near the scale 
k and, evaluated at tree level, describes all quantum effects originating from the high-momentum 
modes. In particular, essentially interpolates between the bare action Sbare — urn fc->oo Tk and 
the standard effective action T = Iihi^qT^. (For more details on the reconstruction of the bare 
action from T^, see [31] .) The fc-dependence of is governed by a FRGE which takes the form 



Here, t = log(k/ko), $ represents all the fields in the theory, and STr denotes a generalized 
functional trace carrying a minus sign for fermionic fields and a factor 2 for complex fields. 
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Lastly, the presence of the cutoff IZk in the FRGE has the remarkable consequence of rendering 
the contribution of STr finite and peaked around p 2 ~ k 2 ; hence, an additional UV regularization 
of the trace is dispensable [TO ] [TTT [31] . 

As the FRGE (TT]) cannot be solved exactly we have to resort to an approximation scheme. A 
standard one is to take an ansatz (truncation) for r&[<l>] of the form 

r k [$] = J2ui(k)om, (2) 

i 

with a finite number of operators 0j[3>] constructed from the fields <£, and with U{ denoting scale- 
dependent coupling constants of mass dimension di. The non-perturbative (3- functions governing 
the RG flow of the Ui's are obtained by substituting ([2]) into (JXJ) and computing the coefficients 
of the operators Oi[$>] appearing on the RHS. The reliability of the results so obtained is checked 
by systematically enlarging the truncation space. 
In the following, we take r^.[$] to be of the form 

r fc [$] =rf + 5 gf + ^ h , (3) 

where T| r , S" gf , and S &h are the gravitational part, the classical background-gauge fixing term, 
and the corresponding ghost action, respectively. The derivative expansion of T| r up to fourth 
order is given by 

rf = J d^g 

where R denotes the Ricci scalar, C 2 = CavpcC^ 1 " 7 is the square of the Weyl tensor, and 
E = R 2 — AR^vRV + R^ U p a RV ypa is the integrand of the topological Gauss-Bonnet term. The 
couplings Ui are related to the dimensionful Newtons constant G and the cosmological constant 
A by uq = A/(8nG), u\ = —1/(16ttG). Truncations including the polynomials of the Ricci scalar 
have already been extensively analyzed in \21 \ I22| , [23l 124] . The inclusion of the tensor structures 
like C 2 , responsible for the key features of higher-derivative gravity, is novel, and provides an 
important test for the reliability of the approximations done so far in the literature. In particular, 
as explained in [32], the inclusion of these terms is necessary in order to show that asymptotic 
safety is not spoiled by those same terms that make gravity perturbatively non-renormalizable. 

The gauge fixing is carried out via the background field method. In the presence of four- 
derivative operators, it is natural to consider a gauge-fixing term which also contains four deriva- 
tives. Following [6] we use 

S gi [g;g} = ^ J d 4 x^F,Y^F u (5) 

with Fp = D u hp V — jDph, Y^ u = D 2 g pv . Here, = g^p — g^ denotes the fluctuations of the 
metric g^ u around the background g^ u , and the background-covariant derivative. Owing to 
the derivatives in Y pLI the ghost sector of the theory contains, in addition to the usual complex 
ghost fields, a "third" ghost [33]. In the sequel, we will work in the limit a — > oo, which leads to 
considerable simplifications. 

The construction of the non-perturbative /3-functions for the couplings ui then proceeds by 
substituting the ansatz ([3|) into ([!]) and projecting the RHS of the equation onto the curvature 



Uq + U\R 



10 n 1 n 6 

— — i? + — C 2 + —E 



3A 



2A 
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monomials contained in (|4|). Going beyond the one-loop approximation [6] thereby requires a 
class of background metrics which is generic enough to disentangle the coefficients multiplying 
R 2 and the tensorial terms, and, most importantly, simple enough to avoid the appearance of non- 
minimal higher-derivative differential operators inside the trace. While the maximally symmetric 
backgrounds used up to now are insufficient in the former respect, a generic compact Einstein 
background £ (which we take without Killing or conformal Killing vectors and without boundary 
for simplicity), satisfying = ^g^u, is sufficient to meet both criteria and allows one to 
determine the non-perturbative /^-functions of the linear combinations 



U 2 



co 9 
3A + 6A ' 



U3 



2A + A 



(6) 



Finding a NGFP for these linear combinations implies that, barring miraculous cancellations, 
none of the couplings oj, 9 goes to infinity, and A and either uj or 9 must be non-zero. 
Substituting ([3]) and equating g = g afterwards, the left-hand side of (pQ) gives 



d t Tf = J d 4 x^[ii + uiR + (ii 2 - gW 3 ) R 2 + u 3 R 



r>fj,upal 
fiupa 1 1 



To evaluate the RHS of (H|) we recast eq. (jU as 



J d 4 x^[u + Ul R + (ua - |n 3 )i? 2 + 2u 3 R pu R' M/ + (n 3 + {)E\ 



(7) 



(8) 



Noting that the variation of E does not contribute to <5 2 r| r /(5<I>5$, the second variation can be 
obtained using the general results [33J , which considerably simplify once the background Einstein 
metric is substituted. Performing the transverse-traceless decomposition of the metric fluctuations 
and ghost fields [M] brings the operator inside the trace into (almost) diagonal form, with the off- 
diagonal terms vanishing in the limit a — > oo. Following \22 \ 123 } 124"]. the Jacobians arising from 
this change of variables can be represented as Gaussian integrals of appropriate auxiliary fields, 
in a similar fashion to the ghost representation of the Faddeev-Popov determinant. The main 
virtue of this decomposition in combination with the choice of background is that all covariant 
derivatives and curvature tensors organize themselves in terms of Lichnerowicz Laplacians, 



A 



-D 



2R«/4> aP , 



-2? 



R 



A 0L <f> = —D 



(9) 



which then commute with all other quantities appearing under the traces. Subsequently, the 
cutoff operators IZ^ are constructed in such a way that the modified propagators are obtained by 
replacing z — ► P n ,k(z) = z + lZ n ^{z) for z = A n £. The resulting flow equation for the truncated 
Tk then becomes 

dtTkl®} = S 2T + Shh + Sit + S , 



(10) 



where 



S2T 



^Tr 2T 



^Tr 



dt{<iu 2 {Pl k -&l L ) + ( : u 1 -2u 2 R)R. Qtk }' 

1 2" 

6u 2 Pg k + {u 1 -2u 2 R)P 0l k + 3 u o 



Sit 
S = 



|Tr 1T 



"|Tr 



dtRi.k 
Pi,k 

dtRo, k 
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Here, u\, = 2u% — U3/3 and the indices "2T", "IT" and "0" indicate that the trace is over 
transverse-traceless matrices, transverse vectors and scalars, respectively. These traces are then 
evaluated using general heat kernel methods adapted to the Lichnerowicz Laplacians. The re- 
sulting expressions are projected onto (H]) evaluated on an Einstein background, giving rise to the 
/?- functions for the dimensionful couplings n.;. As a non-trivial check, we have verified explicitly 
that the expansion of the resulting /^-functions to leading order in A precisely agrees with the 
(projected) universal one-loop result [5J, [6] • 

In order to study the RG properties of the theory it is useful to switch to the dimensionless 
couplings gi = k~ di Ui and consider the /^-functions dtg% = $i- Our particular interest is on 
fixed points of the /3-functions, Pi(gi) = 0, which could provide a non-perturbative definition of 
quantum gravity within the AS program. We find that the /^-functions for the couplings contained 
in (J3J indeed give rise to such a NGFP0 

gl = 0.00442 , g\ = -0.0101 , g\ = 0.00754 , gl = -0.0050 . (12) 

The "universal", i.e., gauge-independent [5], combination GA takes the fixed point value (GA)* = 
0.427, which, together with the value of g%, is in good agreement with previous computations, 
notably the Einstein-Hilbert truncation \17 \ \18 \ fl"9j I20j. the i? 2 -truncation [21], and the trunca- 
tions using a polynomial expansion of /(-R)-gravity [22l [23l [240 • Notice that the finite values for 
g\ and g\ imply a finite value of A*, to be compared with the one- loop result A* = 0. Thus the 
non-perturbative corrections captured by the FRGE shift the fixed point underlying the asymp- 
totic freedom obtained within perturbation theory to the NGFP (|12D featuring in the asymptotic 
safety program. 

An important characteristics of the NGFP are its stability properties. The linearized RG flow 
around the NGFP, dtg% = T$ij(gj — gj), is governed by the stability matrix By = djPi\^ with 
stability coefficients (negative eigenvalues of B) 



O = 2.51, 0i = 1.69, 
and associated normalized eigenvectors 

V = { 0.12 , 0.10 , -0.06 , 0.99} T , 
V 2 ={ 0.74 , -0.46 , 0.48 , -0.11 } T , 



02 = 8.40, 3 = -2.11, (13) 

Vi = { -0.20 , 0.74 , -0.10 , 0.63} T , 
y 3 = { 0.07 , -0.21 , 0.97 , -0.09 } T . 



We observe that the inclusion of the C 2 coupling leads to real stability coefficients at the NGFP. 
This is in contrast to the complex stability coefficients and the corresponding spiraling approach 
of the RG flow characteristic for f(R)-type truncations, and rather reflects the behavior found 
within the perturbative one-loop computation [6] . 

1 For conciseness we only give the numerical results for the optimized cutoff function Rk(z) = (k 2 — z)0(k' 2 — z) 
[35] , but we have verified that the use of a smooth exponential cutoff confirms the picture reported below. 

2 Note that in previous works the background was assumed to be maximally symmetric, which means an Einstein 
space with zero Weyl tensor. In such case R^paR^ 9 " = R 2 /6 and from © it follows that one looses track of M3. 
The value for ^3 constitutes a new result of asymptotically safe quantum gravity. 



5 



Note that the transition from asymptotic freedom to asymptotic safety lifts the degeneracy 
of the marginal couplings. In particular, the negative stability coefficient #3 indicates that the 
corresponding eigendirection is now UV-repulsive. The condition that a RG trajectory is asymp- 
totically safe, i.e., that it approaches the fixed point as k — > 00, implies a relation between the 
couplings Q{. In the neighborhood of the NGFP, this relation allows to express (73 in terms of the 
other couplings, 

g 3 = -0.116 + 0.745#o - 2.44l0i + 11.06g 2 • (15) 

It is remarkable that within our truncation the critical surface remains three-dimensional, as in 
the f(R) case [22l [231 OS] ■ 

A profound consequence of the UV limit of gravity being controlled by a NGFP is the possible 
removal of the poltergeist haunting higher derivative gravity, as proposed in (2J El [27]. In this 
course, we observe that, besides the massless graviton and a massive scalar typical of f{R) theo- 
ries, classical actions of the form (JJJ) give rise to a massive spin-2 poltergeist which is problematic 
for the unitarity of the resulting quantum theory (see, e.g., [HE]). Its mass is given by (minus) 
the pole of the (Euclidean) propagator, m| = A/(167rG). In the quantum regime associated with 
the NGFP, the couplings ui depend on the RG scale, u\ = k 2 g^, 112 = g\ 7^ 0, u 3 = g\^ 0. Thus, 
crucially, A* 7= 0, and m 2 is not constant, but grows as k 2 . This running naturally affects the 
position of the poles in the propagator. Given that Tfc provides an effective description for the 
physics at the scale p 2 ~ k 2 , the mass has to be evaluated at this corresponding momentum scale. 
Identifying p 2 = k 2 , the poltergeist propagator develops a pole at 

k 2 = -m 2 (k 2 ). (16) 

There are then two possibilites. If A* > 0, this effective mass diverges as k — > 00, so that this 
unphysical mode decouples from the theory in the UV. If A* < 0, then, provided the couplings 
enter the linear quantum regime where m 2 {k 2 ) ~ 02k 2 with C2 = —gfX*, before the pole is met, 
the latter will be moved further by the RG running, and indeed completely removed if C2 < — 1, 
viz. Fig. [TJ Unfortunately, our truncation is not precise enough to determine A* and therefore 
draw any definite conclusions. Note, however, that the general mechanism here does not rely 
on our truncation details, but is a generic consequence of the NGFP. The fact the theory is 
asymptotically safe instead of asymptotically free is vital for this argument. 

To summarize, the /3-functions of our higher-derivative truncation possess a non-trivial fixed 
point with three UV-stable and one UV-repulsive eigendirections, in accordance with the asymp- 
totic safety scenario. While a proof of asymptotic safety for gravity would, most likely, include 
the analysis of an infinite series of curvature terms, which is beyond current FRG-methods, the 
mutually consistent results obtained by systematically extending the truncation subspaces un- 
der consideration lends strong support to the idea that gravity is indeed asymptotically safe. 
Moreover, the transition from asymptotic freedom to asymptotic safety provides an interesting 
new perspective on the "poltergeists" appearing in higher-derivative gravity theories, so that 
asymptotically safe quantum gravity may actually be safe from such hauntings. 
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Figure 1: RG-scale dependence of the Poltergeist-mass. In the scaling regime of the NGFP (grey 
region) m 2 , oc k 2 . For the proportionality factor C2 < — 1 the pole of the effective propagator is 
removed dynamically due to the running of the mass. 
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